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Quantum networks hold the promise for revolu- 
tionary advances in information processing with 
quantum resources distributed over remote loca- 
tions via quantum-repeater architectures. Quan- 
tum networks are composed of nodes for storing 
and processing quantum states, and of channels 
for transmitting states between them. The scala- 
bility of such networks relies critically on the abil- 
ity to perform conditional operations on states 
stored in separated quantum memories. Here 
we report the first implementation of such con- 
ditional control of two atomic memories, located 
in distinct apparatuses, which results in a 28-fold 
increase of the probability of simultaneously ob- 
taining a pair of single photons, relative to the 
case without conditional control. As a first ap- 
plication, we demonstrate a high degree of indis- 
tinguishability for remotely generated single pho- 
tons by the observation of destructive interfer- 
ence of their wavepackets. Our results demon- 
strate experimentally a basic principle for en- 
abling scalable quantum networks, with applica- 
tions as well to linear optics quantum computa- 
tion. 

In recent years, it has been established that quantum 
memory is an essential component for the distribution of 
entanglement over arbitrarily long distances using quan- 
tum repeaters [HQ. In the quantum repeater protocol, 
entanglement is distributed by entanglement swapping 
through a chain of spatially-separated entangled pairs of 
particles. Without memory, all pairs need to be entan- 
gled at the same time for the entanglement distribution 
to succeed, an event whose probability decreases expo- 
nentially as the length of the chain increases. On the 
other hand, if it is possible to store the entanglement 
in spatially-separated quantum memories and if one has 
a trigger that unambiguously heralds the entanglement 
once it is achieved 0], it is possible to build up the chain 
of entangled pairs by entangling different parts of the 
chain at different times. By conditioning the evolution 
of the whole system to the output of its different parts, an 
exponential enhancement is attained in the probability of 
success of the protocol, which leads, for example, to the 
possibility of scalable long-distance quantum communi- 
cation JJ- Let us note also that Linear Optical Quan- 
tum Computing [H, [H, 0] or quantum state engineering, 
with schemes working in an iterative manner Q , are both 
damped exponentially by low success probability due to 
the present lack of synchronized single photon sources. 
Real-time control for synchronization of many sources is 



then a promising way to boost the probability of simul- 
taneous generation of many target quantum states, and 
thereby to enable the practical realization of elaborate 
procedures. 

Since photons are the basic carriers of information over 
long distances, the necessity of using memory for com- 
munication implies the need to control the exchange of 
quantum information between matter (stationary qubits) 
and light (flying qubits). Great progress has been 
achieved recently in this direction for different systems 
that could work as single nodes of a distributed quan- 
tum network tflf. For example, generation of photon 
pairs [13, EH, ! storage of single photons [H, [3] , and 
high-efficiency retrieval of stored single excitations [l5[ 
were implemented with atomic ensembles. Entangle- 
ment between atoms and spontaneously emitted light was 
demonstrated with both single trapped atoms [16|, [TtJ 
and atomic ensembles [HI, [l9j]. Heralded entanglement 
between two cold atomic ensembles by means of a sin- 
gle stored excitation was demonstrated last year Q, fol- 
lowed more recently by the achievement of probabilis- 
tic (a posteriori) entanglement between two excitations 
stored in different ensembles 20] . In the continuous vari- 
ables regime, deterministic entanglement has been ob- 
tained between two vapor cells at room temperature [2l| . 

A key point that has not been experimentally ad- 
dressed up to now, however, is the extent to which these 
different systems and techniques enable scalable quan- 
tum networks. Only recently three groups reported the 
use of feedback to enhance the probability of generating a 
photon using a heralded single photon source with mem- 
ory [l9|, H3, HH . In the present work, we report the first 
implementation of real-time conditional control of two 
distant quantum memories. The memory nodes consist 
of ensembles of cold atoms that can each store a sin- 
gle collective excitation in a probabilistic, but heralded, 
way • Since this excitation can be transferred with high 
efficiency to a light field in the single-photon regime [l5| , 
the system functions as a heralded single-photon source, 
ideal for applications to quantum repeaters. The condi- 
tional control allows us to store an excitation in one en- 
semble, while waiting for a trigger signal indicating the 
presence of an excitation in the other ensemble. Relative 
to operation without this conditional control, we attain 
a factor of 28 increase in the probability to generate si- 
multaneously a single photon from each system. Further 
increase is currently limited by the finite coherence time 
of our memory [ill, [24[ • 

For applications in quantum repeaters, it is crucial 
that single photons emitted by two different ensembles 
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FIG. 1: Overview of the experiment, a, Setup for simultaneous storage of single collective atomic excitations in two 
ensembles (L and R). A write pulse initially couples the g — ► e transition, resulting with small probability in the emission of a 
photon (field 1) on the e — ► s transition. Fields 1l,1a are then coupled to polarization-maintaining (PM) fibers and directed to 
detectors D\l and Dir, respectively. In a given trial, if neither of these detectors registers an event, the system is taken back to 
its original condition by means of a strong read pulse coupling the s — > e transition. On the other hand, if one of the detectors 
registers an event, the control electronics uses the Mach-Zehnder intensity modulators (I.M.) in both write and read (panel 
b) pathways to turn off the write-read sequence for the corresponding ensemble. Later, if the other detector also fires, the 
control electronics generates a ready signal, indicating the simultaneous storage of two excitations, b, Setup for simultaneous 
generation of two single photons. The ready signal activates the read intensity modulators, releasing simultaneously the read 
pulses for both ensembles. As a result, two photons (in fields 2l,2r) are generated, and then coupled to the two inputs of 
a PM fiber beam splitter, whose outputs are directed to detectors D2 a ,D2b- A rotatable half-wave plate (A/2) controls the 
polarization of field 2l- 



are indistinguishable when combined at a beam split- 
ter, since this is at the heart of the technique to achieve 
measurement-induced entanglement 0, 0, HIJ and entan- 
glement swapping [4j between remote atomic systems. 
As a first application of our control system, we per- 
form a time-resolved two-photon interference measure- 
ment to quantify the indistinguishability of the generated 
single photons, obtained simultaneously from atomic ex- 
citations stored for different amounts of time. We com- 
bine the photons at a beam splitter and observe a non- 
classical suppression in the rate of joint detections at 
its two output ports (2fjj . This suppression is the re- 
sult of a destructive two-photon interference, as first 
demonstrated in a parametric down conversion system 
by Hong, Ou, and Mandel [27[. More recently such sup- 
pression was observed for photons emitted successively 
by a single source 0, [2i| or simultaneously by different 
sources, such as spatially separated down-converters [30| 



and trapped atoms [3l|. Our results show a suppres- 
sion of (77 ±6)% for the probability of having two pho- 
tons leave the beam splitter through different ports, from 
which we infer an overlap £ ~ 0.90 for the wavepackets 
of the two single photons. 

Our experimental setup is sketched in Fig. [1] (see also 
Methods section). The two ensembles consist of pencil- 
shaped clouds of cold cesium atoms located in two dif- 
ferent vacuum chambers, 2.7 m apart. In the beginning 
of each trial, all atoms are optically pumped to the hy- 
perfine ground state \g) (see Fig. QJ,). The atomic level 
structure for the writing process consists of the initial 
ground state \g) (6Si/ 2 ,F = 4 level of atomic cesium), 
the ground state |s) for storing a collective spin flip 
(651/2, F = 3), and the excited level |e) (6P 3/2 ,F = 4). 
Then a weak write pulse, lasting 38 ns, excites the g — > e 
transition. With a small probability qi ~ 0.005 <C 1, 
the atomic ensemble spontaneously emits a photon (field 
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1) on the e — > s transition, into the solid angle of our 
detection system. For our experimental conditions, the 
detection of this first photon heralds the storage of an 
excitation in a collective, symmetric mode of the whole 
ensemble 0, HEI13- 

This collective excitation can then be retrieved with 
high probability [J, [l5| by a strong read pulse (38 ns 
long and resonant to the s — » e transition) counterprop- 
agating with respect to the write beam; see Fig.[]i>. The 
read pulse results in the generation of a second photon 
(field 2) in the direction opposite to field 1 [ll[ . For both 
ensembles, a detection in field 1 occurs with probability 
Pi = 0.12% and is followed by a detection in field 2 with 
conditional probability p c ~ 8.5%, which corresponds, 
after taking the losses in the field-2 channels into ac- 
count, to q c ~ 34% retrieval efficiency for the collective 
mode at the output of the ensemble [HI]. If no detec- 
tion in field 1 is registered in a given trial, the read pulse 
is fired in order to optically pump the atoms back to 
their initial condition. For the chosen pi, the normalized 
intensity cross-correlation function between fields 1 and 
2 is measured to be approximately 512 — 23 for both en- 
sembles [HI], corresponding to a field 2 well within 
the single-photon regime, with a large suppression of its 
two-photon component [15j. For our system, 512 > 2 
is already a strong indication of nonclassical correlations 
between fields 1 and 2 [l(|. The parameter w quantifies 
this suppression by examining the probability of gener- 
ating two photons in the same pulse, normalized by an 
equally bright Poisson-distributed source [Hj]. Classical 
fields must satisfy the Cauchy-Schwarz inequality w > 1; 
for independent coherent states, w = 1, while for thermal 
fields, w = 2 [331 ]. The parametric dependence of w on 
gi2 for our system was investigated in detail in Refs. [TH 
and [H, from which we infer w ~ 0.17 for 512 = 23. 

Fields 1_l,1_r are guided to the respective detectors 
{D 1L ,D 1R ) (Fig. QJ,). The ficld-2 outputs of the two 
ensembles, on the other hand, are combined at a fiber 
beam splitter, whose outputs are then directed to two 
detectors D^a and Z?2b (Fig. [TJs). To exploit the quan- 
tum memory to speed up probabilistic quantum protocols 
that require concurrent state preparation in two atomic 
ensembles, we have designed and implemented a custom 
logic circuit that allows conditional control of the writ- 
ing, storing, and reading operations for the atomic exci- 
tations. Upon receipt of a field- 1 detection signal from 
one ensemble, the circuit gates off the write and read 
pulses for that ensemble , thereby storing a collective ex- 
citation in the atoms [la . Il9l . l33j ] . The write-read pulse 
train in the other ensemble is not affected. The storage 
stops and the excitations in the ensembles are read out 
when either of the following events occurs: (1) A field- 1 
photon is detected in the second system, prompting the 
circuit to release read pulses into the first and second 
ensembles, thereby simultaneously retrieving the stored 
excitations from both ensembles; (2) A pre-determincd 
maximal storage time Ai max set in the circuit is reached 
in the first ensemble. Then, the storage ceases, and the 



ensembles are optically pumped back to the original con- 
dition. The logic circuit returns to its dormant state, 
passing all the write and read pulses to the ensembles, 
until the next field- 1 detection signal triggers its func- 
tion. 




FIG. 2: Probabilities pu and P1122 of coincidence detection 
as functions of the number N of trials waited between the 
independent preparations of the two ensembles (L, R) with 1 
excitation each. Filled squares give the joint probability pu 
of simultaneously preparing the two ensembles. Open circles 
give the joint probability P1122 of preparing the two ensem- 
bles and detecting a pair of photons, one in each output of 
the beam splitter, in fields 2l,2r. Error bars indicate ±VC 
photon counting noise, where C is the number of counts. To- 
tal number of counts for the evaluation of pu was 234244, 
and 614 for pn22- Total number of trials was 3.36 x 10 9 , cor- 
responding to 6 hours of data taking. The polarizations for 
fields 2l,2r were set to be orthogonal. 

Our logic circuit is designed specifically to increase pu, 
the probability of having excitations stored simultane- 
ously in each ensemble. Given an initial event in one 
ensemble, the filled squares in Fig. [5] show how pu in- 
creases as a function of the number of trials N of duration 
Atrial — 525 ns that we wait for the other ensemble. The 
values in this curve were obtained by performing the ex- 
periment using a maximum number of trials N max = 23 
(corresponding to our specific Ai max = 12ps), recording 
a file with the whole history of events, counting the num- 
ber of events where the second ensemble was prepared up 
to N trials after the first one, and then dividing this num- 
ber by the total number of trials N t . After 23 pulses, we 
observe an increase T\\ — 44 inpn, close to the expected 
value of T\\ = (2N max — 1) = 45 for the case pi << 1 
(see Methods Section). Note that for our experimental 
conditions, the conditional logic reduces (with reference 
to a try-until-success strategy) the number of trials N t in 
a negligible way, by up to 6% (i.e., 2p\N max — 0.055). It 
is also important to point out that the advantage of such 
conditional logic will be (exponentially) more pronounced 
for larger systems involving more than two memories. 

Our interest, however, is not directly in pu, but rather 
in the probability P1122 of obtaining joint detections for 
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field 2 after the two systems (L, R) have each stored a 
single excitation. In Fig. OP1122 (open circles) is given 
as a function of N, for the case where the fields 2^,2^ 
from the two ensembles have orthogonal polarizations. 
Evidently, P1122 does not increase by the same factor as 
does pn, but grows instead only up to J-1122 = 28 [i.e., 
28 times the value obtained without the control circuit 
(N = 0)]. This behavior is expected from the fact that 
the retrieval efficiency for the stored excitation in the 
ensemble that is first prepared decays until the other en- 
semble is prepared and the read pulses are fired. The 
coherence time r c for our system is about 10 [is [H^HH, 
and can also be directly inferred from the decay with N 
of the probabilities for one and two photo-detections in 
field 2 (see Appendix). 

To characterize the indistinguishability of the photons 
in fields 2^ and 2^, we perform a measurement of the sup- 
pression of joint-detection events at detectors £>2a, -D2& in 
the same trial. In Fig. [3J we observe that the conditional 
joint-detection probability P22 decreases for the situation 
where fields 2^,2^ are combined with the same polariza- 
tions (filled squares), compared to the case of orthogo- 
nal polarizations (open circles). The detection times td 
are obtained from the record of events in our acquisi- 
tion card, which has 2 ns resolution, much smaller than 
the duration of the photon wavepackets. From this list 
of detection times, we obtain the time difference t be- 
tween the two detections. To quantify the coincidence 
suppression shown in Fig. [3J we consider the visibility 
V = (p± — P\\)/p±, where p± gives P22 with orthogonal 
polarization and p || gives that with parallel polarization. 
Note that V — 1 for ideal single photons with perfect 
overlap of their wavepackets, while V = for completely 
distinguishable fields. Obtaining p± and p || from the in- 
tegration over r of the respective data points in Fig. [3J 
we find V = 0.77 ± 0.06 (see also Appendix). 

For photons orthogonally polarized, we fit our mea- 
surements for the joint-detection probability P2 = 2 ( r ) m 
Fig. El to a Gaussian (dashed line), 
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We find T = 18.4 ± 0.2 ns, from which the photon du- 
ration can be inferred assuming identical wavepackets. 
For perfectly transform-limited photons, the suppressed 
P2 2 (r) with parallel polarizations can be obtained from 
the one with perpendicular polarizations by multiplica- 
tion by a constant scaling factor / = 1 — Vfu- When one 
introduces frequency jitter Aw, P^k 7 ) can be expressed 
aslil [Hill: 



[1 - V fit cos(A^r)] (2) 
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The green line gives such a fit with Vfit — 0.80 ± 0.02 
and Aui/2tt = 4±4MHz. The jitter error bar is obtained 
by doubling the x 2 fitting parameter. This simple fit 
agrees well with the data and leads to a time-bandwidth 
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FIG. 3: Conditional joint-detection probability P22C 7 ") of 
recording events in both D2 a and D26, once the two en- 
sembles are ready to fire, as a function of the time differ- 
ence t between the two detections. Filled squares (empty 
circles) provide the results when the field-2 outputs of the 
two ensembles are combined with parallel (orthogonal) po- 
larizations. The dashed line is obtained from a Gaussian fit 
of the orthogonal-polarization data, with half- width at 1/e 
equal to T = 18.4 ± 0.2 ns. The solid line is obtained from 
the dashed one by multiplying it by 1 — Vfu cos(Aa>r), with 
Vfu = 0.80 and = AM Hz. Error bars indicate ±VC 

photon counting noise, where C is the number of counts. 



product equal to 1.2 ± 0.2, providing a good indication 
that our photons are close to transform limited. Oth- 
ers measurements are necessary to investigate this issue 
further. 

The main cause of visibility reduction V < 1 is that the 
two-photon components for the conditional fields 2^,2^ 
for both ensembles are necessarily nonzero (lo| . From the 
inferred value of two-photon suppression w = 0.17, we 
estimate with a simple model that the maximal achiev- 
able visibility for perfect overlap between the two fields 
is V„iax = 0.85 (see Appendix). By comparing our mea- 
sured visibility V = 0.77 to V maX7 we infer that the over- 
lap of the field-2 wavepackets is £ ~ 0.90, where £ = 1 
for perfect mode-matching. The overlap mismatch can 
be partly explained by a non-zero polarization extinc- 
tion ratio in the polarization-maintaining fibers (-14 dB). 
The visibility also decreases due to a small misalign- 
ment introduced by the rotation of the half-wave plate 
to switch between the two polarizations. The alignment 
was checked by analyzing the events where the two pho- 
tons are created in different readout intervals, thus ex- 
hibiting no interference (see Appendix). In this case, 
the orthogonal-polarization joint-detection level is 0.08 
smaller than for parallel polarization, which should in- 
crease the visibility by 0.02. Another small decrease in 
the visibility should come from the measured imbalance 
(0.51/0.49) of the fiber beam splitter. Note that the vis- 
ibility of the two-photon interference can in principle be 
increased by reducing the intensity of the write pulses 
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(and thereby q\ and hence w jl5l.l33|). but this comes at 
the expense of a reduced counting rate. 
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FIG. 4: Probability densities pd for the wavepackets of fields 
2l and 2r derived from the distribution of detection events 
over time td- Filled squares (empty circles) give the condi- 
tional wavepacket for field 2l (2jj), for detections occurring 
only in the same trial as a field-I detection. The solid (dashed) 
line gives the unconditional wavepacket for ensemble L (R). 
The data-taking duration was 200 s for ensemble L, and 300 s 
for ensemble R. All curves are normalized by their areas. Er- 
ror bars indicate statistical errors. 

We also performed independent measurements to as- 
sess the temporal overlap of the two wavepackets. The 
temporal shape of the conditional wavepackets of fields 
2l,2r are given, respectively, by the filled squares and 
empty circles in Fig.HJ These wavepackets were obtained 
by blocking the fields from the other ensemble, disabling 
the logic circuit, and considering only the field 2 detec- 
tions occuring in the same trial as a field 1 detection. The 
two wavepackets are quite similar, consistent with the 
good suppression we measured, and both have temporal 
widths around T c ~ 13 ns (Gaussian fit), which is also 
consistent with the expected width T/v2 = 13 ± 0.2ns 
extracted from Fig. El by considering identical wavepack- 
ets. In Fig. 2] we also show the unconditional wavepackets 
for fields 2^,2^. They are considerably different from the 
conditional ones, indicating that part of the light we col- 
lect from the ensembles is not correlated to field 1, but 
likely arises from processes not related to the collective 
atomic state [Hf. 

In summary, we have demonstrated a 28-fold enhance- 
ment in the probability for simultaneous generation of 
single photons from two remote atomic ensembles, as a 
result of the real-time control of the storage of collective 
excitations in the ensembles. We have used this enhance- 
ment to demonstrate the indistinguishability of two pho- 
tons generated by independent systems, through the ob- 
servation of a destructive interference of their wavepack- 
ets that results in the coalescence of the two photons 
at a beam splitter. From this measurement, we have 



also inferred that the generated single photons are nar- 
rowband, and their wavepackets are close to transform- 
limited. We emphasize that absent conditional control 
of the two remote systems, or without quantum memory, 
data of comparable quality to that presented here would 
have required continuous acquisition over more than two 
weeks, which is prohibitive. The fundamental control 
for quantum state manipulation implemented here will 
be integral to future advances with networks of quantum 
memories, including for quantum repeaters. Our work 
thereby paves the way to scalable quantum networks over 
distances much longer than set by fiber optic attenuation. 

METHODS 

Experimental details. Magneto-optical traps are used 
to form the clouds of atoms, and are switched off for 6 ms 
every 25 ms period. After waiting for the trap magnetic 
field to decay 24| , a train of write and read pulses ex- 
cite the sample during the last 2 ms. The write pulse 
is 10 MHz red detuned from the g — > e transition. The 
transverse waist of the write beam is 200 /xm, and its 
peak power P wr ite ~ 2/xW. We collect the light emit- 
ted by the ensemble in a polarization-maintaining (PM) 
fiber, whose projected mode on the ensemble corresponds 
to a beam with 50 fim. waist intersecting the write-beam 
direction at a three-degree angle jl5j |. In the experiment, 
the read pulse is delayed from the write pulse by 300 
ns, leaving time for the pulses to be gated off after the 
heralding signal, which occurs 100 ns after the write pulse 
due to propagation delays. 

Increase in probability. Assume that p\ gives the 
probability per trial of storing a collective excitation, and 
that it is possible to wait up to N trials before reading 
out the excitation and releasing the corresponding single 
photon. The probability of having two ensembles storing 
excitations in the same trial is then 

Pn = Pi[pi + 2 [(1 ~Pi)pi + (1 -Pifpi 

+ --- + (l-pi) N - 1 Pi]} 
w (2N - l)p\ , when p 1 « 1 . 

The factor of two in the above expression accounts for 
the two possible orders in which the ensembles can be 
prepared. 
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APPENDIX A: DECOHERENCE 

Figure [5] shows the variation of the conditional prob- 
abilities of detecting one (p|) and two (p£ 2 ) photons in 
fields 2^ and 2#, once the two ensembles are ready to fire, 
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as functions of the number TV of trials that occur between 
the two field-1 detections. This figure was obtained from 
the same raw data as Fig. [21 Fields 2^ and 2^ have then 
orthogonal polarizations, and are combined at a beam 
splitter as shown in Fig.[Tj3. In order to obtain the quan- 
tities in Fig. [5l we divided the number of coincidences in 
field 2 that followed two temporally separated detections 
in field 1 by the number of times the two ensembles were 
prepared with that specific time separation. The solid 
lines are fittings considering an exponential decay of the 
conditional probability p c for the second photon from ei- 
ther of the two ensembles, once a detection has occurred 
in field 1. We assumed the same p c and decay time for 
the two systems. Note that the two systems were actually 
set up to have similar p c and similar Raman linewidths 
for transitions between the hyperfine ground states [24|, 
which should correspond to the system coherence time. 
The expressions used to fit were then 

P22(N) = P le- N ' N ° , (Al) 

pi(N) = P±±P^ p§2{N) . (A2) 

We assume above that p£ 2 always involves one photon 
coming from an excitation stored during N trials. In 
this way, we are neglecting the two-photon component 
of field 2, as well as diverse sources of background. For 
P2 , the first two terms take into account that the con- 
ditioned detection event can be originated from either 
an ensemble that has just been excited, or a stored (for 
N trials) excitation. The third term subtracts the prob- 
ability of having a joint detection on field 2 (p 2 c gives 
the probability of detecting an event on one detector 
and zero on the other). From the fitting, we then ob- 
tain p c — 0.091 and N c — 18, corresponding to a co- 
herence time r c = N c x 0.525/js = 9.5 /is. Keeping in 
mind the simplicity of the above expressions, which do 
not take into account any background in field 2 or its two- 
photon component, the inferred conditional probability 
p c is then consistent with the independently measured 
values of about 0.085 for each ensemble. 



p n22 (N) = pi{pip 2 c 2 (0) + 2 [(1 - Pl ) PlP ^(l) + (1 - 



A plot of this expression, considering the p c and iV c ob- 
tained from the fits in Fig. [SJ is shown as the red curve in 
Fig. [6] The quite reasonable agreement with the experi- 
mental data (filled squares) indicates then that the exper- 
imentally observed increase in P1122 can be understood 
by the increase in pn provided by the circuit combined 



From the above discussion, it is then straightforward to 
obtain an expression taking into account decoherence for 
the measured P1122, presented in Fig. [2] Note first that, 
in the ideal case of very long coherence time (N c —►00), 
P1122 can be obtained from the p\\ expression presented 
in the Methods section by multipling it by the conditional 
probability p 22 ^) °f obtaining a pair of photons with 
effectively zero delay (N/N c — > 0) between them: 




5 10 15 20 

N 

FIG. 5: Conditional probabilities p% and P22 of measuring one 
(red circles and green triangles) and two (black squares) pho- 
tons in field 2, respectively, once the two ensembles are ready 
to fire. The red and black curves are fits using Eqs. (|AT|) 
and l|A2|) . as discussed in the text. The green line is 0.95 x 
the red line, as the field 2 level measured by D26 is always 5% 
lower than for D2 a , indicating a possible difference of detec- 
tion efficiency of this magnitude. 



Considering the measured value of p\ and the value of 
Pc obtained from the fits of Fig. [5j we then obtain the 
green line plotted in Fig. [6l In order to introduce deco- 
herence in this analysis, each term of the p\\ expression 
deduced in the Methods section should be multiplied by 
the proper p 22 (N) as defined in Eq. lAlt 



) W 2 C 2 (2) + • • • + (1 - Pif-'piP^N - 1)] } . (A4) 
I 

with the decoherence of the stored collective excitation. 

Note finally that P1122 times the number of trials per 
second gives the rate of conditional joint detections in 
fields 2 L , 2 R . In this way, from Fig. [5] we can see that a 
larger coherence time could still enhance this detection 
rate by up to a factor of 1.6 for N — 23 (enhance the 
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factor ^1122 from 28 to 45). An increase on the condi- 
tional probability p c would also greatly improve this rate, 
since it scales with p 2 c . As discussed in the text, we infer 
that the probability q c of extracting the photon from the 
ensemble is about 34% for our experimental conditions. 
Thus, considering the same amount of losses on the field 
2 pathways and the same detection efficiencies, we infer 
that ideally, if one achieved q c — 1, p c can still be in- 
creased by up to a factor of 3, which would increase the 
joint detections rate by 9. This indicates a good prospect 
for further optimizations of our system in the future. 
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FIG. 6: Probability P1122 of coincidence detection as functions 
of the number N of trials waited between the independent 
preparations of the two ensembles (L, R) with 1 excitation 
each. Filled squares give the measured joint probability P1122 
of preparing the two ensembles and detecting a pair of pho- 
tons, one in each output of the beam splitter, in fields 2l,2r. 
Error bars indicate statistical errors. The polarizations for 
fields 2l,2_r were set to be orthogonal. These experimental 
results were also presented in Fig. [2] The green curve gives 
the theoretically expected Pnll' f° r the ideal case of very long 
coherence time, as given by Eq. IA3I The red curve gives the 
theoretical P1122 for the case of a finite coherence time, Eq, IA4l 



APPENDIX B: VISIBILITY AS A FUNCTION OF 



In the case of two indistinguishable single-photon 
wavepackets combined at a 50/50 beam splitter (BS), 
no coincidences should be observed at the two ouputs of 
the BS . However, if the combined fields are not per- 
fectly overlapping single photons, coincident detections 
at the output of the BS can occur due to the two-photon 
component in each input port. In this section, we evalu- 
ate the loss of visibility due to this effect with a simple 
model. 

Suppose that the two ensembles are prepared each with 
stored excitation, as heralded by a detection in field 1 for 
both ensembles. Let's denote P n the probability of find- 
ing n photons in field 2, and assume, for simplicity, the 



various P n are the same for both ensembles. In each field 
(before the BS), the two-photon suppression is charac- 
terized by the parameter w [33] : 



w 



PI 



(Bl) 



so that the two-photon component can be written as: 

wPl 



Pi 



(B2) 



Let us now combine the two fields at the BS. The proba- 
bility to have one photon at each output of the BS, when 
the two wavepackets do not overlap (e.g., if they have 
orthogonal polarizations) is given by: 



P-L = 



P2 
2 



P2 
2 



wP{ wP{ 



El 

2 



(B3) 



where the two first terms corresponds to the terms with 
two photons in one input mode of the BS, and the third 
term to the case with one photon in each input mode. 
The factor 1/2 corresponds to the 50% chance that the 
photons split at the BS. In this simplified calculation, we 
neglect the case where we have two photons in one input 
and one in the other one, whose probability is on the 
order of Pf. 

If the two fields overlap perfectly at the BS (parallel 
polarizations), the term with one photon in each input 
does not lead to coincidences, and the probability to have 
one photon in each output is then: 



P\\ 



wPf 



wPf 



(B4) 



Taking Eqs. (|B3[) and (|B4|) into account, we find that the 
visibility can be written as: 



V 



P-L ~P\\ 
P± 



1 



f 



(B5) 



In our case, we have 1712 ~ 23 for the two ensembles, 
from which we estimate w w 0.17 [HI]. This leads to a 
maximal visibility of V max = 0.85 for a perfect overlap 
£ = 1.0 between the fields. From our measured visibility 
of 0.77, we then estimate an overlap £ = 0.90. 



APPENDIX C: JOINT-DETECTION LEVELS 
FOR EVENTS IN DIFFERENT TRIALS 

In Fig. [TId, we show how the conditional probability 
of detecting two photons, when the (L, R) systems are 
ready, decreases as a function of the delay between the 
two detections for the situation where the fields 2^,2^ 
are combined with the same (red) or orthogonal (black) 
polarizations. It corresponds then to Fig. [3] The time t& 
of the detections is obtained from the recording of events 
in our acquisition card, and it refers to a fixed reference 
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that marks the beginning of the 525 ns repetition periods. 
From this list of detection times, we obtain the relative 
delay r when the two detections occur in the same trial. 




-40 40 -40 40 -40 40 

x (ns) x (ns) x (ns) 



viously done in section II, let's denote Pi and P2 the 
probabilities of finding respectively one or two photons 
in each field 2. For the sake of simplicity, these proba- 
bilities are taken equal for the two ensembles. Including 
two-photon events, the probability for coincidence in the 
two detectors is given for the center peak (Fig. 2b) by: 

\pI + P2 (CI) 

The first term takes into account the two cases where 
the single photons are both reflected or transmitted at 
the beam splitter. The second one corresponds to the 
case where two photons arriving at one input of the beam 
splitter are split into the two arms. Higher-order cases, 
which involve for instance two photons in each input, 
or two photon in one input and one in the other, are 
neglected. 

1.0 1 1 



FIG. 7: Conditional joint-detection probability and D2b, once 
a ready signal is generated as a result of the two ensembles 
being ready to fire, as a function of the time difference r 
between the two detections, (b) The two detections occur 
within the same trial, (a) Detector D^a fires first and D2b 
fires after the next ready signal, (c) Same as (a), but with 
the detector order inverted. The red circles (black squares) 
provide the results for field 2 from the two ensembles having 
orthogonal (parallel) polarizations. 

We also show in Fig. the cases where the two detec- 
tions in Di a and Du, occur in different trials, with the 
event in one detector occuring when the two ensembles 
are ready and the event in the other detector occuring 
the next time the ensembles are ready. Figures [7^, and[7Ji 
give the cases in which one detector or the other regis- 
ters an event first. The fact that the signal level is similar 
in both cases, with different polarizations, indicates that 
there is no large misalignment when the half-wave plate 
is turned to switch between the two polarization con- 
figurations. Even though, if we integrate the curves in 
Figs. [Th. and [7b, the value obtained for orthogonal po- 
larization is about 0.08 lower than the one obtained for 
parallel polarization. We confirmed this value by cal- 
culating it also from other different-trials peaks (for de- 
tection events separated by up to 5 ready signals). The 
curves with different polarizations were taking at alter- 
nate cycles of half-hour data taking, just turning a single 
half- wave plate between them. In this way, we believe the 
decrease for orthogonal polarizations comes just from a 
small misalignment in the fiber input introduced by this 
operation. 

The level of the peak in (b) obtained with orthogonal 
polarizations should be half that observed in (a),(c) for 
the case where pure single photons arrive at the beam 
splitter. The experimental observed ratio r is found to 
be r — 0.60 ± 0.05. This ratio can be explained by the 
two-photon component of our generated state. As pre- 




x ( ns) 

w 

FIG. 8: Visibility V of the two-photon coincidence suppres- 
sion as a function of the integration window for the time dif- 
ference between the two field-2 detections. The integration is 
made from — t w to t w . 



In a similar way, the probability for the others peaks, 
where detections occur in different trials, can be written 
as 



{p 1 + 2P 2 y 



(C2) 



This expression corresponds basically to the mean photon 
number going to one detector times the mean photon 
number going to the other one. 

With P2 = wPi/2, and by neglecting higher-order 
terms, the ratio becomes: 



1 I 



2 1 + 2wPi 



(C3) 



With Pi = 0.085 and w — 0.17, this expression gives then 
an expected value r = 57%, which is consistent with the 
observed one. 
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APPENDIX D: VISIBILITY AS A FUNCTION OF 
INTEGRATION WINDOWS 

In Fig. [8] we show the results of the measurement of V 
for different integration windows around r = 0. We see 
that for an integration window around the center, from 
t = —6 ns to r — 6 ns, the visibility is 80±10%, while the 
integration using the whole window gives V — 77 ± 6%, 
indicating that the suppression is roughly uniform for all 
t, which is also consistent with having close to transform- 
limited wavepackets for both fields 2 L ,2 R . 



sitioned around the center of the respective wavepackets. 
In order to analyze the conditional field- 2 wavepackets 
overlap in Fig. [3J and also in the above Figs. [7| and [4j we 
introduced in the analysis an additional time window, 
only 44 ns long around the conditional field 2. As can 
be seen in Fig. |4j this corresponds to consider the whole 
conditional ficld-2 wavepackets. 



APPENDIX E: TIME WINDOWS 



The electronic time windows used for field 1 and field 
2 detections were 80 ns and 90 ns long, respectively, po- 
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